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Charge glass in two-dimensional arrays of capacitively coupled grains with random 

offset charges. 
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V) ■ We study the effect of random offset charges in the insulator to conductor transition in systems of 

capacitively coupled grains, as realized in two-dimensional arrays of ultrasmall Josephson junctions. 
In presence of disorder, the conductive transition and charge ordering at nonzero gate voltages are 
both destroyed for any degree of disorder a at finite temperatures T, in the thermodynamic limit, 
but crossover effects will dominate at length scales smaller than f CTj T = T~ v f(o-T~ v ) , where v is 
the thermal critical exponent of the zero-temperature charge glass transition. The conductance is 

. — ' linear and thermally activated but nonlinear behavior sets in at a crossover voltage which decreases 

as temperature decreases. For large disorder, the results are supported by Monte Carlo dynamics 
* ■ simulations of a Coulomb gas with offset charges and are consistent with the thermally activated 



behavior found in recent experiments. 

ently prepared systems, and therefore different degrees 
of disorder, can still lead to the same activated behavior 
of the linear conductance is still lacking. In particular, 
measurable nonlinear effects resulting from disorder have 
not been considered so far. 

In this work, we consider the effects of random off- 
set charges on the charge-anticharge unbinding transi- 
tion in an otherwise perfect two dimensional array of ca- 
pacitively coupled normal or superconducting grains in 
the classical regime of large charging energy E c and fi- 
nite temperature T. We show that, in the presence of 
offset charges, the unbinding transition and charge or- 
dering at nonzero gate voltages are both destroyed for 
any degree of disorder a in the thermodynamic limit, 
consistent with the common behavior of different exper- 
imental arrays. The conductance is linear and thermally 
activated but crossover effects will dominate at length 
scales smaller than a disorder and temperature depen- 
dent length, £,{j : T = T~ u f(aT~ 1 '), where v is the thermal 
critical exponent characterizing the T — charge glass 
transition. Nonlinear behavior sets in at a characteristic 
voltage V c ~ T/^ a ,T- For larger disorder, the results are 
supported by a Monte Carlo simulation of the nonlinear 
conductance leading to V c ~ T 1+u with a very rough 
estimate v ~ 1.7. 

The electrostatic energy of the net mobile charges Qi = 
e*n,j, multiples of an elementary charge e*, located at the 
sites i of an array is given byfcrO 

E=l^2(Q i + q t )Cr j \Q t+qt ) (1) 

where e* = e, the electron charge, for a normal grain 
and e* = 2e for a superconducting grain. C$j is the ca- 
pacitance matrix and represents a net offset charge 
in a grain induced by charged impurities trapped near 
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Two-dimensional arrays of Josephson junctions where 
the charging energy of the grains is much larger than 
the Josephson coupling energy are interesting systems 
where collective charging-eflEacts can be studied exper- 
imentally in great detailtmoQ. In these systems, the 
net charges in the grains have a long range logarith- 
I ' mic interaction when the junction capacitance C is small 
. enough. The charges play the same role as vortices in 
' the resistive behavior of the array, when capacitive ef- 
fects can be neglected. This allows for the possibility of 
a Kosterlitz- Thouless (KT) charge unbinding transition 
corresponding to an insulating to conducting transition 
at finite temperature! Below the critical temperature 
| T c , neutral charge pairs would be present in the insulat- 
ing phase while above T c , screening of the electrostatic 
potential gives rise to free isolated charges and a conduct- 
ing phase. Experimentally, a transition has in fact been 
observed in the conductance of small oaaaptance arrays 
of normal and superconducting grainsa-BB at an appar- 
ent temperature consistent with the estimate based on 
logarithmically interacting charges. Recently, however, 
a closer analysis of the experimental data of different 
groups, revealed that the onset of finite conductance at 
finite temperatures could be described just as well by 
thermal activation of free chargesBEl with an activation 
energy E a ss \E C when the grains are in the normal state 
and E a jE c + A for superconducting grains, where 2A 
is the superconducting energy gap at zero temperature. 
This interpretation neglects the logarithmic-iateraction 
between charges expected for an ideal arrayM and sug- 
gests in turn that, in the experimental systems, the in- 
teraction is screened by some unknown mechanism and is 
essentially short ranged. Although disorder effects, in the 
form of random offset charges, which may be trapped near 
the grain-substrate interfacetl, could be a possible expla- 
nation for these results, an understanding of how differ- 
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the grain-substrate interfaceoO. The equilibrium proper- 
ties in the classical limit are described by the partition 
function Z = J2{n } e~ E ^ kT , where the sum is over all 
integers rii and the offset charges qi act as quenched dis- 
order with a nonzero uncorrelated component. We also 
assume that tunneling of charges between the grains pro- 
vides the dynamics to equilibrate the charge distribution 
and the junction dissipation is sufficiently smallEJ. Since 
an offset charge equal to an integer number of e* , can be 
compensated by shifting rij — > — m, one can consider 
offset charges with a probability distribution of width 
a < e* peaked at zero. In general, the elements of the 
inverse capacitance matrix C^ 1 fall off exponentially for 
Tij > A, where A ~ ^/CjC g is a screening length, C 
the capacitance between neighboring grains and C g is 
the capacitance to the ground. We consider the regime 
r/A << I where! C*^ 1 = -^logry - 1/4C. Under 
these conditions, the interaction energy above describes 
a two dimensional Coulomb gas of logarithmically inter- 
acting charges Qi in a background of quenched random 
charges qi. This problem has been considered before in 
relation to the effects of geometrical disorder in classical 
superconducting arrays in a magnetic-field where Qi rep- 
resents vortices and qi a random fluxa and in the gauge- 
glass modefl of disordered superconductors where the 
are random but correlated. 

In the absence of offset charges, qi = 0, a charge- 
anticharge unbinding transition in the KT universal- 
ity class has been predicted!!! at a finite temperature 
kT c ps e* 2 /8nC. However, in the presence of offset 
charges with a nonzero component of uncorrelated dis- 
order, < qiqj > = a 2 Sij, the behavior is drastically 
changed as bound charge- ant icharge pairs are unstable 
at a sufficiently large length scale £ determined by the 
degree of disorder. In fact, the net charge fluctuation 
due to random charges in a domain of area irL 2 will neu- 
tralize a test charge Q — e* when L = = e* /^/ira. 
For length scales larger than charges are unbound 
and will contribute to the linear conductance G provided 
that charges are mobile. This shows that even for small 
a, disorder is a strongly relevant perturbation. 

In general, there will be also a contribution from cor- 
related disorder. However, this is either irrelevant at the 
fixed point of the pure system or at most marginal for suf- 
ficiently small disorder. To lowest order in the wavevec- 
tor k, the moments of disorder distribution are given by 
< q k >= and < q k q-k > = c 2 + erffc 2 + 0(/c 4 ). The 
first term corresponds to uncorrelated randomness as dis- 
cussed above and the second term to offset charges corre- 
lated as random dipoles p = q \n —Tj\ , where r, and Vj 
are neighboring sites, since < PkP~k >— v\k 2 . Higher 
order terms correspond to multipole correlations. Us- 
ing standard methods to convert .-the partition function 
of Eq.(I) to a sine-Gordon modelo and the replica trick 
to average the free energy over disorder, we obtain an 
effective reduced Hamiltonian in terms of n replicas 

H/kT = \ f d 2 r[± £(V^) 2 + (^) 2 £ tftf 



+(^?) 2 E V ^ V ^ -2/E cos ( 2 ^ Q )]' ( 2 ) 

where K — E c /2n 2 kT. Higher order gradient terms have 
been neglected. In the absence of disorder this model 
has a massless phase for irK > 2 when the charge fugac- 
ity y is irrelevant corresponding to an insulating phase. 
For ttK < 2, y is relevant leading to a conducting phase. 
From renormalization- group arguments, the relevance of 
the disorder perturbation on the massless phase can now 
be determined from the eigenvalue A = d — x where 2x 
is correlation function exponent evaluated at the unper- 
turbed fixed point. For uncorrelated disorder, a ^ 0, it 
is sufficient to consider the diagonal contribution a = (3. 
This has an eigenvalue A = 2 in the massless phase since 
x = (d — 2) = in two dimensions, which corresponds to 
a scaling behavior a' 2 = a 2 b 2 under a change of lattice 
spacing by a factor of 6, in agreement with the domain 
argument given above leading to a characteristic disorder 
length £ CT cx 1/a. For correlated disorder alone, cr 2 ^ 0, 
the contribution from the diagonal term ^ a V</? a V<£> Q 
just renormalizes the value of K but the offdiagonal per- 
turbation J2 a ^(3 ^J^a^Jtpp is marginal, A = . The 
higher order gradient terms in the Hamiltonian are then 
irrelevant. 

The disorder dependent length £ CT is the scale up to 
which the charges Qi are mobile. Beyond this scale, they 
are neutralized by the random background of quenched 
charges qi. At a scale L > the system is better de- 
scribed as a set of pinned charges interacting by a short 
range potential and we expect that it behaves similarly 
to a charge glass. Since these charges play the same role 
as vortices in a two-dimensional superconductor, one ex- 
pects that glass order occurs at T = only, with a ther- 
mal correlation length £*r oc T~ v where v is a thermal 
critical exponent associated with the T = fixed pointfl. 
Assuming that the model of Eq.(l) is in the same uni- 
versality class, we expect that v > 2. Standard crossover 
scaling leads to a correlation length = T~ u f(<rT~ v ) 
where the scaling function fix) is a function of the di- 
mensionless ratio x = ^t/C<t with f(x) ~ I for x » x c 
and f(x) ~ 1/x for x << x c so that £ CTj t oc 1/a when 
oT~ v < x c and £ CT:T cx T~ v for oT~ v > x c . For an infi- 
nite system, the predicted insulating -conducting transi- 
tion in the absence of disorder is destroyed by any finite 
density of uncorrelated offset charges but crossover ef- 
fects may lead to an apparent transition if is larger 
than the system size. Note that the above arguments are 
not strictly correct in the limit a = or T = as it does 
not allow for the possibility of a transition to a true in- 
sulating phase of bound pairs ofxharges at finite T when 
(7 = or for quantum tunnelingtl at T = 0. In principle, 
such effects could be incorporated in the crossover scaling 
analysis at the expense, «f some major extra complica- 
tions but experimentsaBB indicate that there is a range 
of temperatures over which purely classical behavior is 
observed and that behavior characteristic of the absence 
of disorder is never observed. The crossover scaling anal- 
ysis outlined above should be relevant to experimental 
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systems at temperatures above the crossover tempera- 
ture where quantum tunneling dominates™. 

In addition to £, another length scale, £v-nis set by an 
applied finite voltage V across the systemE3. The addi- 
tional contribution to the energy has the form J^. QiExi 
for an uniform electric field E in the x direction. Thermal 
fluctuations alone, of typical energy kT, leads to a charac- 
teristic length t;v ~ kT / e*E over which single charge mo- 
tion is possible. For £y > £, when charges are unbound, 
this leads to a linear response to the applied voltage and 
a finite conductance limited by the tunneling rate across 
each junction in the array. For £y < £, the contribution 
from the applied voltage must be balanced against the in- 
teraction energy with other charges leading to nonlinear 
behavior. A crossover from linear to nonlinear behavior 
occurs when these two lengths are comparable, i.e., at 
a characteristic voltage V c ~ T/£ a , T = T u + 1 g{oT~ v ), 
where g{x) — 1/ f(x) . 

The finite correlation length £ in presence of disor- 
der also determines the behavior of the activation en- 
ergy in the regime of linear conductance. For vanishing 
voltage, the current arises from thermal dissociation of 
the most weakly bound charge pairs. For normal met- 
als, the energy necessary to unbind a typical pair is 
E = log(£/a)+E , where E a ~ E c /2 is the core energy 
of a pair with one lattice spacing separation. For suffi- 
ciently large disorder such that £ ~ a, the conductance 
should then have-an Arrhenius behavior with an activa- 
tion energy givenclby E a = E /2. However, for moderate 
disorder the activation energy has a temperature and dis- 
order dependent logarithmic correction log(£ CTl r/a). In 
the pure case where er = and £ a j? - * oo, this implies 
an infinite activation energy and the linear conductance 
vanishes consistent with V c ~ T/£ CT ,T — ¥ 0. The loga- 
rithmic dependence on of the activation energy also 
indicates that its value is only drastically changed for 
small disorder or low temperatures where £ CTi t is large 
but is roughly insensitive to larger values of disorder or 
temperature. For superconducting grains, the arguments 
above still apply except that an additional energycl equal 
to the superconducting energy gap E g is required to cre- 
ate quasiparticles which allow charges on the grains to 
tunnel to neighboring grains, in the absence of Joseph- 
son tunneling. 

In the presence of a gate voltage V g on each grain, in 
the absence of disorder there will be an induced charge 
CgVg on each grain acting as a frustration parameter^/ — 
CgVg/e* analogous to the flux quanta per plaquette in 
Josephson junction arrays in a magnetic field. The prop- 
erties of the system should be periodic in / with period 
1. In particular, the Coulomb gap width oscillates as a 
function of V g with period e*/C g , often used as an esti- 
mate of C g . For fractional values of / the ground state 
is a charged ordered state with an average number / of 
elementary charges per grain. For rational / = p/q, this 
leads to a ground state with a unit cell of size qx q and a 
discrete degenerate ground stateO. A finite temperature 
transition is expected from the charge ordered phase to 
a disordered phase in addition to the charge-anticharge 



unbinding transition. In presence of disorder, the behav- 
ior at integer / will be the same as / = 0, i.e., at any 
finite temperature the interacting charge system is disor- 
dered and the conductance is finite. This behavior should 
persist for fractional values of / since, in addition to the 
discrete symmetry excitations such as domain walls, ther- 
mally excited net charges interacting logarithmically are 
also present which will be screened by disorder at length 
scales larger than £ leading to a linear thermally acti- 
vated conductance. However, one still may question if 
charge order remains stable in presence of disorder. It is 
sufficient to consider the case £-^j= 1/2 since the ground 
states for / < 1/2 are less stableEil. At f — 1/2 and in the 
absence of disorder, it is known that the ground stateMl 
consists of 2 x 2 unit cells with an antiferromagnetic ar- 
rangement of n, — / = ±1/2 net fractional charges, in 
units of the elementary charge e* . The lowest energy 
excitations are domain walls separating the two degener- 
ate ground states, ±1/4 charges at domain wall corners 
and unit excess charges interacting logarithmically. In a 
coarse-grained description, the variables in a unit cell can 
be replaced by Ising variables S r = ±1, which describe 
the charge ordered state, -and charges Qr = SV/4 av- 
eraged over rthe unit cel£± This leads to an effective 
Hamiltonianl 

H/kT = L x ^ S r s r> + Zk 2 K Qr G(R - R') Q R , 

<r,r/> R.R> 

+L 2 Y / S r q r + 4TT 2 K^2Q R G(R-r)q r (3) 

r rR 

where G(r) = log(r/a) + ir/2. The Ising variables S r 
are coupled antiferromagnetically by a nearest neighbor 
interaction L\ and to a random field proportional to q r . 
The excess charges Qr are still coupled to the random 
offset charges q r and will therefore lead to finite conduc- 
tance as discussed before. Since the offset charges cou- 
pled as random fields to the Ising variables, charge order 
will be destroyed even for small disorcLer. There will be 
however a characteristic length scaleE3 (f a ~ exp(c/u 2 ), 
beyond which random offset charges destroy the charge 
ordered state. For sufficiently small a, £ CT > £ CT and there 
can be an intermediate length scale where charge order 
remains but excess charge is effectively unbound. Similar 
results should apply for f < 1/2 and we expect that dis- 
order destroys charge ordering for any finite gate voltage. 

We have studied the linear and nonlinear conductance 
of the Coulomb gas in Eq.(l) numerically using classical 
Monte Carlo nonequilibrium dynamicsQ, regarding the 
Monte Carlo time as real time. This should be a reason- 
able approximation in the overdamped limit and when 
quantum effects can be ignored. It also assumes that the 
tunneling rate which provides the-equilibrium charge dis- 
tribution satisfies the global ruleEJ when it depends on 
the energy difference of the whole array. In the simu- 
lations, a Monte Carlo step consists of adding a dipole 
of unit charges and unit length to a nearest neighbor 
charge pair, (n^rij), corresponding to a tunneling of the 
positive or negative charge by a unit of length. The ex- 
ternal electric field biases the added dipole, leading to 
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a current as the net flow of charge is in the direction 
of the electric field. The unit of time corresponds to a 
complete Monte Carlo pass through the system. The 
linear conductance Gl = Hthe^qJ/E is obtained from 
the equilibrium current fluctuations, without imposing 
a voltage bias, from the fluctuation-dissipation relation 

G L = wr! dt < I W( t ) >■ 
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FIG. 1. (a) Linear conductance Gl with zero gate voltage 
(/ = 0) and (b) with nonzero gate voltage (/ = 1/2). The 
straight lines indicate an Arrhenius behavior. Temperature T 
is measured units of E c /2tv 2 . (c) Nonlinear conductance J/E 
as a function of temperature, for a — 0.5 and (/ = 0). (d) 
Scaling plot of the nonlinear conductance in (c). 

Figs. ( la, b) show the behavior of the linear conductance 
Gl in the absence of disorder, a = 0, and in the pres- 
ence of a (large) random Gaussian distribution of offset 
charges with a standard deviation a = 0.5. Since large 
disorder leads to small £, we can use a small system size 
L x L with L = 16. When a = 0, the conductance ap- 
pears to vanish at T c w 1.4 for / = and sa 0.8 for 
/ = 1/2, corresponding to the critical temperature for 
the conductive transition for zero and nonzero gate volt- 
ages. However, for a = 0.5 , the conductance remains 
finite much below the critical temperature for the pure 
system. In fact, as indicated in the figure, the data in 
the low temperature range can be reasonably fitted to 
the Arrhenius activated form G = G e~ Ea ^ kT in agree- 
ment with a disordered phase. Moreover, the calculation 
of Gl from current fluctuations ensures that the acti- 
vated behavior is not an artifact of a voltage bias since 
it is obtained without an external electric field. In ex- 
periments, however, a voltage bias is often used and non- 
linear effects should be taken into account. In Fig.(lc) 
we show the behavior of the nonlinear conductance as 
a function of the applied electric field and temperature 
for a = 0.5 and / = 0. For small E, J/E tends to 
the constant finite value of Fig. (la) which decreases with 
temperature. For the highest temperature, T = 0.8, the 



range of E in which this behavior is apparent is more 
pronounced. For increasing E it clearly crosses over to a 
nonlinear behavior in agreement with the analysis lead- 
ing to a crossover voltage V c which is both temperature 
and disorder dependent. The nonlinear conductance can 
be cast into a scaling fornxl, J/EGl = F(E/E C ), where 
F is some scaling function with F(0) — 1. For large dis- 
order we have E c ~ T v+X , and a scaling plot allows an 
estimate of v . This is shown in Fig. (Id) where a rea- 
sonable data collapse is obtained by adjusting E a and v 
for the largest temperatures, giving v ~ 1.7, in rough 
agreement with the T = stiffness exponent for gauge 
glass models!! This implies that the characteristic volt- 
age V c where nonlinearities set in should exhibit a non- 
trivial power law V c ~ t u+1 and should be accessible 
experimentally. 
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